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One challenge of using millimeter-wave radar for measuring clouds is achieving 
system calibration. The delicate nature of millimeter-wave components and the harsh 
environment in which they operate may cause undetected changes in the system response 
unless regular system calibration is performed. Scattering properties of the ocean surface 
have been widely used as a calibration reference for airborne and spaceborne microwave 
sensors. However, at millimeter-wave frequencies, the ocean surface backscattering 
mechanism is still not well understood, in part, due to the lack of experimental 
measurements. During the Cirrus Regional Study of Tropical Anvils and Cirrus Layers- 
Florida Area Cirrus Experiment (CRYSTAL-FACE), measurements of ocean surface 
backscattering were made using a 94-GHz (W-band) cloud radar onboard a NASA ER-2 
high- altitude aircraft. The measurement set includes the normalized ocean surface cross 
section over a range of the incidence angles under a variety of wind conditions. Analysis 
of the radar measurements shows good agreement with a quasi-specular scattering model. 
This unprecedented dataset enhances our knowledge about the ocean surface scattering 
mechanism at 94 GHz. The results of this work support the proposition of using the ocean 
surface as a calibration reference for airborne millimeter-wave cloud radars and for the 
ongoing NASA CloudSat mission, which will use a 94-GHz spaceborne cloud radar for 
global cloud measurements. 
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ABSTRACT 


Scattering properties of the ocean surface have been widely used as a calibration reference 
for airborne and spacebome microwave sensors. However, at. millimeter- wave frequencies, the 
ocean surface backscattering mechanism is still not well understood, in part, due to the lack of 
experimental measurements. During the Cirrus Regional Study of Tropical Anvils and Cirrus 
Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE), measurements of ocean surface 
backscattering were made using a 94-GHz (W-band) cloud radar onboard a NASA ER-2 high- 
altitude aircraft. The measurement set includes the normalized ocean surface cross section over a 
range of the incidence angles under a variety of wind conditions. Analysis of the radar 
measurements shows good agreement with a quasi-specular scattering model. This unprecedented 
dataset enhances our knowledge about the ocean surface scattering mechanism at 94 GHz. The 
results of this work support the proposition of using the ocean surface as a calibration reference for 
airborne millimeter-wave cloud radars and for the ongoing NASA CloudSat mission, which will 
use a 94-GHz spacebome cloud radar for global cloud measurements. 
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1. Introduction 


Clouds piay a critical role in the Earth's climate system. The vertical structure and spatial 
distributions of clouds are important in determining the Earth's radiation budgets, which affect 
global circulations and ultimately climate. However, the lack of fine-scale cloud data is apparent in 
current climate model simulations (Houghton et al. 1995; Stephens et al. 1990). In cloud detection, 
millimeter-wave cloud radars have gained favor because of their high scattering efficiency, low 
power consumption, and compact size. A number of airborne millimeter-wave cloud radars have 
been developed (Pazmany et al. 1994; Sadowy et al. 1997; Li et al. 2003). Meanwhile, a 94-GHz 
spacebome cloud radar is in preparation for the NASA CloudSat program (Stephens et al. 2002). 
Another proposed spacebome W-band cloud radar is the joint Japanese/European Earth Cloud, 
Aerosol and Radiation Explorer (EarthCARE) mission radar, as described by the European Space 
Agency (2001). 

One challenge of using millimeter-wave radar for measuring clouds is achieving system 
calibration. The delicate nature of millimeter-wave components and the harsh environment in 
which they operate may cause undetected changes in the system response unless regular system 
calibration is performed. The calibration uncertainty specification of the CloudSat 94 GHz cloud 
radar is 2 dB (with a goal of 1.5 dB) (Stephens et al. 2002). Typically, a radar may utilize internal 
circuitry to monitor the variation of the power levels of the radio frequency (RF) transmitter and the 
drift of system gain, but characterization of the antenna, front-end waveguides, and the interfaces 
between these components, are not included with internal calibration. External calibrations using 
point targets, such as a trihedral comer reflector, have been performed for ground-based radars (Li 
et al. 2003; Sekelsky 2002). However, it is difficult to calibrate airborne or spacebome radars using 
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point targets due to the difficulty in separating the return of the calibration reference from 
background clutter. 

The ocean surface has been widely used as a calibration target for airborne and spacebome 
microwave radars and radiometers. Numerous studies of ocean surface scattering have been 
performed at microwave frequencies (Jones etal. 1976; Plant 1977; Valenzuela 1978; Masuko et al. 
1986). Based on experiences with operating millimeter-wave cloud radars on the ground and on 
airborne platforms, it will be essential to perform CloudSat radar calibration checks periodically 
once it is launched. The ocean surface will be a valuable calibration reference since there will be 
many measurement opportunities under clear-sky conditions. However, in part because of the lack 
of experimental measurements at millimeter-wave frequencies, the ocean surface backscattering 
mechanism is still not well understood. Meanwhile, attenuation caused by water-vapor and oxygen 
absorption in the lower troposphere is significant at millimeter- wave frequencies. It is therefore 
necessary to correct this attenuation to reduce the uncertainty of the calibration. In addition to the 
radar calibration, the ocean surface return can also be used in estimating the path- integrated 
attenuation (PIA) along the radar beam. The PIA then is used as a constraint for retrieving radar 
attenuation, rain rate, or cloud particle microphysical properties (Meneghini et al. 1983; Iguchi and 
Meneghini 1994). 

During July 2002, the Cloud Radar System (CRS), a 94-GHz (W-band) pulsed-Doppler 
polarimetrc radar developed by NASA Goddard Space Flight Center, was operated on a NASA ER- 
2 high-altitude (nominal 20 km) research aircraft in support of the NASA CRYSTAL-FACE 
program (Jensen et al. 2003). During the same experiment, the 9.6-GHz ER-2 Doppler Radar 
(EDOP) (Heymsfield et al. 1996), dropsondes (Hock and Franklin 1999) as well as other remote 
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sensors were also operated on the ER-2. Table 1 shows CRS system parameters during CRYSTAL- 
FACE. The radar calibration was performed using a trihedral comer reflector after the CRYSTAL- 
FACE . The calibration result was verified by side-by-side inter-comparison with the University of 
Massachusetts (UMass) ground-based Cloud Profiling Radar System (CPRS) 95-GHz cloud radar 
(Sekelsky and McIntosh 1996), which has been well maintained and calibrated over the past 
decade. Inter-comparison measurements between the CRS and CPRS for similar cloud volumes 
were conducted during autumn 2002. Figure 1 shows comparison of CRS and CPRS W-band radar 
reflectivity profiles. The mean difference of the profiles in clouds is within 1 dB. The details of the 
CRS hardware, installation on ER-2, system calibration, and preliminary cloud measurements from 
CRYSTAL-FACE are described by Li et al. (2003) and McGill et al. (2003). During the 
CRYSTAL-FACE, cean surface scattering was measured by CRS over the southern Florida 
offshore region. The radar measurements are used to investigate the efficacy of using ocean surface 
returns as a reference for calibrating airborne and spacebome W-band radars. 

This paper presents a comparative analysis of the CRS CRYSTAL-FACE data with an 
ocean surface scattering model using three parameterizations of the mean surface slope. In Section 
2 the ocean surface backscatter model is introduced, and the dependency of the ocean surface radar 
cross section on incidence angle and surface wind speed is discussed. The need to correct for 
atmospheric attenuation at 94 GHz is discussed as well. Section 3 explains configuration of the 
instrumentation during the CRYSTAL-FACE and describes the dataset used in this study. Section 4 
presents CRS measurements and compares these measurements with EDOP data and the 
backscatters models. 

2. Background 
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a. Ocean surface backscatter model 


In general, the normalized ocean surface scattering cross section o a is a function of radar 
wavelength, radar beam incidence angle, polarization, ocean surface wind speed, and wind 
direction. For incidence angle 6 smaller than 15°, cr 0 is dominated by large-scale surface waves, and 
at microwave frequencies the quasi-specular scattering theory has been shown to work well in this 
region (Valenzuela 1978; Barrick 1974; Brown 1990). When 6 is larger than 15°, Bragg scattering 
produced by small-scale waves becomes more significant, and therefore, two-scale scatter models 
have been used since they take into account both quasi-specular scattering and Bragg scattering 
(Brown 1978; Chan and Fung 1977). For the CRS operating on the NASA ER-2 and for the 
CloudSat radar, the primary measurement objective is to obtain vertical profiles of cloud and 
precipitation layers. Therefore, these radars are designed to operate at low incidence angles and the 
quasi-specular scattering theory is considered valid. Under this theory, the ocean surface is 
assumed isotropic and the surface wave distribution probability density is only a function of the 
surface mean-square slope, s(v); o 0 is then approximately given as (Valenzuela 1978; Brown 
1990): 


<J o (0,v,X) 


=J 

s(v) 2 cos 4 (0) 


tan 2 (ff) . 
rfv) 2 h 


( 1 ) 


where X is the radar wavelength, v is surface wind speed in ms' 1 , and s(v) 2 is the effective mean- 
square surface slope. The ocean surface effective Fresnel reflection coefficient at normal incidence, 
r e (0,X) = C e (n(X )- 1 .0 )/(n(X )+ 1 .0), and n(X) is the complex refractive index for seawater. For the 
CRS operating frequency, n(X = 3mm) = 3.36-/1.93 at 20 C (Meneghini and Kozu 1990). The 
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Fresnel reflection coefficient correction factor, C e , will be discussed in Section 3c. It is worth 
noting that this quasi-specular can not resolve the dependency of cr 0 on wind direction. 

The effective mean-square surface slope, s(vf, is given by different empirical relationships. 
Cox and Munk (1954) developed a linear relationship based on classical optical scattering data. 
They showed that s(vf = 0.003+5.08xl0' 3 v, varied linearly with wind speed for both a low-pass 
filtered surface where only gravity waves are present, and a “clean” surface which includes 
capillary wave scale roughness. Wu (1972, 1990) reanalyzed the data obtained by Cox and Munk 
(1954) and showed that the mean-square slope varies approximately with the logarithm of wind 
speed. Therefore, s(vf is expressed as s( v f - wo+wj/ugiofv), where wo and w\ are empirically 
determined constants, but different in two wind speed regimes: for v < 7.0 ms' 1 , w 0 = 0.009 and wj = 
0.0276; for 7.0 <v < 20.0 ms' 1 , wo = -0.084 and wj = 0.138. Based on the statistical analysis of the 
a 0 measurements obtained by the Tropical Rainfall Measuring Mission (TRMM) Ku-band (13.8- 
GHz) precipitation radar and surface wind measurements obtained by the passive TRMM 
Microwave Imager, Frelich and Vanhoff (2003) derived wo = 0.0036 and wi = 0.028 for 1.0 <v < 
10.0 ms' 1 , and w 0 = -0.0184 and wi = 0.05 for 10.0 <v < 20.0 ms' 1 . In the following, these three 
surface slope relationships are discussed and compared with radar measurements from CRYSTAL- 
FACE. 


Figure 2a shows er 0 at 94 GHz versus incidence angle at 2.0, 6.0 and 10.0 ms' 1 wind speed 
using the above three different s(vf empirical relationships. In general, as the surface wind speed 
increases, the surface roughness increases and the shape of the <r 0 versus incidence angle curve 
becomes flatter. At a given wind speed, cr 0 calculated using TRMM s(v) 2 decreases more rapidly 
with the increase of incidence angle. For 10.0 ms 1 wind, the results from Cox and Munk (1954), 
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and Wu (1972, 1990) are basically identical, so that these two curves are indistinguishable in Fig. 
2a. Although the a 0 versus 0 relationships are different for the three different empirical models, the 
difference tends to be the smallest near a 10° incidence angle. 

Figure 2b gives a 0 versus wind speed at 5°, 10°, and 15° incidence angles for different s(v) 2 
empirical relationships. For 5.0 <v < 18.0 ms' 1 , a Q calculated using Cox and Munk’s (1954) and 
Wu’s (1972, 1990) s(vf relationships are nearly the same at all three incidence angles except for a 
1 dB discrepancy near 7.0 ms 1 at 5° and 15° incidence angles. This is due to the discontinuity in 
Wu s relationship at a 7.0 ms 1 wind speed. It is worth noting that er 0 is less sensitive to wind speed 
near a 10° incidence angle for all three s(vf relationships when wind speed is higher than 2 ms' 1 . 
This has been reported by many authors from previous measurements using microwave radars 
(Jones et al. 1976; Plant 1977; Masuko et al. 1986). In addition. Figure 2b also indicates that at a 
10° incidence angle, a 0 estimates using the three different s(vf empirical relationships agree within 
1 dB for 3.0 <v < 10.0 ms' 1 . This fact is very useful for calibration of airborne or spacebome radars 
since pointing the radar beam at incidence angles near 10° over an ocean surface minimizes 
uncertainties due to surface wind variations. 

b. Radar equation of ocean surface scattering 

The basic form of a radar equation for surface scatter is given by (Kozu 1995): 

P - 

r SUMnXljJh 2 ’ () 

where, P r = power at the receiver (mW), 
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P t = peak transmit power (kW), 

G a = antenna gain, 

A = radar signal wavelength (m), 

fl, <p = antenna 3 dB beamwidth in horizontal and vertical (radians). 

It = loss between the antenna and receiver port, 
l a = loss between the transmitter and the antenna port, 
latnf= one-way path-integrated atmospheric attenuation 
6 = radar beam incidence angle (radians), 

°o = normalized ocean surface radar cross section, and 
h = altitude of the aircraft (m). 

PG ^ 

The product of a only depends on the radar system parameters. This product is related 

.r^tx 

to the radar constant such as that used by Sekelsky (2002), 

r - 10241n2/t 2 / tt Z r 10 21 

P t G'CJT*Tj3<p ’ ^ 

where c is the speed of light (3xl0 8 ms' 1 ), and r is the RF pulse width (s). 

R c can be evaluated from: 1) measurements of individual parameters in (3); 2) an external 
calibration using a target with known radar cross section, such as a trihedral comer reflector 
(Sekelsky 2002); and 3) the ocean surface. For our case, R c for the CRS was obtained from a series 
of external calibrations using a trihedral comer reflector (Li et al. 2003). Assuming that a beam- 
filled condition is satisfied, cr 0 can be calculated as 
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a. 


(4) 


Px^Lj-rRh 2 

r aim c 

2i 4 10 21 cos(0) ' 

On the other hand, if the ocean surface conditions and cr 0 are known, the radar constant can 
be derived from (4) so that radar calibration is achieved. In Section 4, observational estimates of a 0 
and their dependence on the radar beam incidence angle and surface wind conditions are examined. 

c. Attenuation due to atmospheric absorption 

Water vapor and oxygen molecules do not produce significant backscatter at 94 GHz. 
However, water- vapor and oxygen absorption at millimeter-wave frequencies is much stronger than 
at microwave frequencies (Lhennitte 1988; Clothiaux et al. 1995; Li et al. 2001). Because water 
vapor and oxygen are highly concentrated in the lower troposphere, ocean surface measurements 
made from an airborne or spacebome radar are attenuated by their presence. For tropical and mid- 
latitude regions the attenuation can be significant. In subtropical ocean areas such as Florida, the 
two-way path-integrated attenuation from water-vapor and oxygen absorption can be as high as 7 
dB as a result of the high humidity and high temperature. A practical atmospheric millimeter-wave 
propagation model developed by Liebe (1985) predicts attenuation and path delay of moist air for 
frequencies up to 1000 GHz. The input variables of this model are height distributions (0-30 km) of 
pressure, temperature, and humidity along the propagation path. Using profiles of relative humidity, 
temperature, and pressure measured by the ER-2 dropsondes and radiosondes launched from the 
surface, the attenuation has been estimated and corrected using liebe’ s (1985) millimeter-wave 
propagation model. 

3. Data description and processing 
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a. Radar configuration during CRYSTAL-FACE 

During CRYSTAL-FACE, the CRS was installed in the tail cone of the ER-2 right wing 
superpod, while EDOP was operated from the nose of the ER-2. The measurement geometry of the 
CRS during the CRYSTAL-FACE experiment is illustrated in Fig. 3. Since the CRS was 
configured in a fixed nadir pointing mode, incidence angle scanning in the cross-track direction 
was achieved when the ER-2 made routine maneuvers. When the ER-2 changed its flight direction, 
it rolled its body to one side and gradually changed its heading direction, then rolled back to a level 
position to complete the turn. The incidence angle of the radar beam, 8, is derived using aircraft 
pitch-and-roll angles. 

The CRS antenna has an elliptical beam with a full-width-half-power beamwidth of 0.6° 
crosstrack and 0.8° alongtrack. The ER-2 flew at a nominal 20-km altitude; therefore the CRS 
footprint size on the surface is 210 m cross track and 280 m along track. The RF pulse width of the 
CRS was 1.0 ps (150 m in range), and data were averaged for 0.5 s before being stored onto a solid- 
state recorder. With the 1.0 gs pulse width, the surface footprint is beam-filled up to a 32° incidence 
angle. During the CRYSTAL-FACE flights, the maximum radar beam incidence angle was less 
than 30° when the aircraft made turns. Therefore, all CRS measurements from this experiment were 
valid for the beam-filled condition. Meanwhile, the CRS data were sampled at a 0.25-ps time 
interval (37.5-m range spacing), and the ocean surface return was oversampled by a factor of four. 
Oversampling the surface return is a technique used for airborne and satellite observations which 
significantly reduces the error of the surface return power (Kozu 2000; Caylor et al. 1997). 
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The incidence angles of the CRS radar beam are determined using aircraft pitch-and-roll 
angles, which are provided by the aircraft navigation system at an 8-Hz rate and simultaneously 
recorded along with the 2-Hz radar data. The accuracy of the aircraft pitch and roll is 0.05°. Ideally, 
the radar beam should be pointed along nadir. However, during installation, alignment errors result 
in offsets of the radar beam from nadir in both pitch and roll. These offsets are estimated and 
corrected using o 0 and Doppler velocity measurements of the ocean surface. The offset in pitch 
(along track) was obtained using the Doppler velocity of the ocean surface, which should be 0 ms" 1 
on average. Using this method, the CRS antenna pointing uncertainty was determined to be about 
0.2° in pitch. The offset in roll (cross track) was estimated by comparing er 0 measurements from 
right and left turns. At low incidence angles, <x 0 versus 6 curves measured from right and left turns 
should agree with each other since a D is not sensitive to wind direction at low incidence angles 
(Jones et al. 1976). The CRS beam offset in roll was determined to be less than 0.4° (see Section 
4a and Fig. 10). 

In this study, measurements of <j 0 made by CRS are compared to the measurements marfe by 
the 9.6-GHz EDOP, which has two beams, one pointing to nadir and the other pointing to 33° 
forward from nadir. The EDOP has been well calibrated using the TRMM precipitation radar and 
ocean surface return (Heymsfield et al. 1996, 2000). The beamwidth of the EDOP nadir beam 
antenna is 2.9°, about four times that of the CRS beamwidth. The RF pulse width and sampled gate 
spacing of EDOP are 0.5 ps (75 m in range) and 37.5 m, respectively. For this operating mode, the 
beam filled approximation is only valid for incidence angles smaller than 5°. Therefore, <r 0 
measured by the EDOP nadir beam is performed partially beam-filled correction (Equation [16] in 
Kozu [1995]) so that it is comparable to the CRS data. In order to reduce the errors in the 
measurement of a 0 , EDOP ocean surface return was oversampled by a factor of two. 
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b. Atmospheric attenuation estimation and surface wind speed 


Sounding data profiles, used to estimate atmospheric attenuation, were obtained by the ER- 
2 GPS dropsondes (RD-93 model) developed by the National Center for Atmospheric Research 
(NCAR) (Hock and Franklin 1999). They provided pressure, temperature, humidity, and surface 
horizontal wind, etc. The measurement accuracies of these dropsondes are: pressure, ±1.0 hPa; 
temperature ±0.2 C; humidity, ±7%; horizontal wind, ±0.1 ms' 1 . During the CRYSTAL-FACE, 
between 4 and 8 dropsondes were typically launched during each 5-hour flight. Figure 4 shows 
temperature, relative humidity, and pressure profiles measured by a dropsonde from clear-sky 
conditions at 1954 UTC, 9 July 2002. The latitude and longitude of the dropsonde launch point 
were N23.83° and W86.15 0 , respectively. The estimated two-way integrated attenuation due to 
water-vapor and oxygen absorption is shown in Fig. 4d. For this case, the maximum two-way 
attenuation due to water-vapor and oxygen absorption is 6.7 dB. 

The ER-2 flew 11 science missions during the experiment in which about 50 dropsondes 
were launched. Figure 5a shows the two-way path-integrated attenuation due to water-vapor and 
oxygen absorption under clear weather condition. The results show that the averaged two-way 
integrated water-vapor and oxygen attenuation is approximately 5.8 dB, while the standard 
deviation is 0.55 dB. The high attenuation and large standard deviation reveal the importance of the 
correction of ocean surface measurements for atmospheric attenuation. In addition to the near- 
surface measurements made by the dropsondes, ocean surface wind measurements were obtained 
from nearby surface buoys as well. Figure 5b shows the surface wind speed measured by ER-2 


13 


dropsondes during CRYSTAL-FACE. The measurements indicate that the ocean surface wind 
speed was between 1 and 7 ms' 1 for most ER-2 flights. 

C. Estimate of correction factor C e 

For a smooth surface, the Fresnel reflection coefficient at normal incidence is given by the 
classic formula as E(0,X) = (n(X)-1.0)/(n(X)+l.Q), where n(X) is the complex refractive index of the 
surface materials (Ulaby et al. 1981). At 94 GHz, n(X = 3mm ) = 3.36-/1.93 for 20 C seawater 
(Meneghini and Kozu 1990); therefore, \T(0, X = 3mm)f = 0.409 for a smooth ocean surface. 
However, the ocean surface is generally roughened by gravity waves, surface winds, and 
precipitation. According to Jackson et al. (1992), small-scale surface roughness diffracts the 
incident radiation. This diffractive process reduces the Fresnel reflection coefficient. It is therefore 
necessary to use effective Fresnel reflection coefficient r e (0, X)- C e (n(X) -1.0)/(n(X)+1.0), where C e 
accounts for the surface roughness effects. C e is generally smaller than 1 and is equal to 1 only for a 
perfectly smooth surface. 

Based on the statistical analysis of the ocean surface backscattering measurements obtained 
by the Tropical Rainfall Measuring Mission Ku-band (13.8 GHz) spacebome radar, Frelich and 
Vanhoff (2003) estimated that C t is approximately 0.89 for surface wind speed between 1.5 ms' 1 
and 15 ms' 1 . At 94 GHz, estimates of C e have been difficult to obtain due to the lack of 
measurements of the ocean surface. Here, we compare o 0 calculated from (1) and <x 0 measured by 
the CRS to estimate C e . At low incidence angle (near nadir), c 0 is not useful for estimating C e 
because it is sensitive to different s(vf empirical relationships (see Fig. 2) and has relatively large 
variability with surface wind speed at nadir. However, at a 10° incidence angle, o 0 is insensitive to 



surface wind speed and to s(v) 2 empirical relationships for 3.0 <v < 10.0 ms' 1 . The mean value of 
<t 0 calculated within this wind speed range (3.0 <v < 10.0 ms' 1 ) for all three s(vf relationships is 
6.94 +201og(C e ) dB. <r 0 measurements at 10° incidence angles obtained by the CRS from a different 
turn events have a mean value of 5.85 dB and a standard deviation of 0.6 dB. Since o 0 calculated 
from the model should agree with the measurement, C e is then estimated as 0.88 with an 
uncertainty of 0.16. Errors in system calibration, o 0 calculation, and atmospheric attenuation 
estimates are possible sources of the uncertainty in C e . This C e estimate is almost the same as the 
TRMM estimate, but it does not necessarily mean that C e is frequency independent. More 
experimental data are necessary to examine the dependency of C e on radar frequencies. 

4. a 0 observed by the CRS 

a. <r 0 versus incidence angle 

On 9 July 2002, the ER-2 flew a tropical cirrus mission from its base at Key West, Honda, 
through the Yucatan Channel, then south-southeast into the Caribbean Sea (Fig. 6). On the return 
trip, the ER-2 made a turn through point B (1935 UTC, N21.72 0 , W86.ll 0 ) under clear weather. 
Figure 7 shows the aircraft navigation data, incidence angle, and surface return power. It indicates 
that the ER-2 rolled its body to the right side during the first half of the turn, and rolled back in the 
second half of the turn to resume a level position. During the turns, the CRS radar beam scanned 
along the cross-track direction away from and back to nadir. Figure 8 shows <r 0 measured by CRS 
versus the incidence angle compared with theoretical curves. The a 0 measured by EDOP is also 
shown in Fig. 8 as diamonds. Figure 8 shows a larger variation in EDOP measurements, which 
were oversampled by a factor of two than the CRS measurements, which was over sampled by a 
factor of four. This observation is consistent with the results reported by Caylor et al. (1997) and 
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Kozu et al. (1995), and indicates that the errors in a 0 estimates can be reduced by oversampling 
along the radar range. 

It is difficult to obtain the exact surface wind conditions at turn point B, but measurements 
were made by dropsondes launched at 1601 UTC (N23.75 0 , W86.16 0 ), and at 1644 UTC (N19.01 0 , 
W86.91 0 ) on the outbound trip, and at 1954 UTC (N23.83°, W86.16 0 ) on the return trip (see 
locations in Fig 6). The near surface wind speeds measured by these dropsondes were 2.6 ms' 1 , 3.1 
ms' 1 , and 6.8 ms' 1 , respectively. Results using the models described in Section 2 with wind speeds 
of 2.5 ms' 1 , 3.0 ms' 1 and 6.5 ms" 1 are plotted for comparison in Fig. 8. In Fig. 8a, the CRS 
measurements fall between the 2.5 ms' 1 and 6.5 ms' 1 Cox and Munk model predictions, and they are 
closest to the 3.0 ms' 1 curve. However, near a 10° incidence angle, the CRS measurements (~ 6 dB) 
agree well with the model results for the three wind speeds. For EDOP, only the 3.0 ms" 1 Cox and 
Munk model result is plotted (dashed-dotted line. Fig. 8a). Figure 8b shows that both the CRS and 
EDOP measurements match well with Wu’s (1990) surface slope relationship for 3.0 ms" 1 wind 
speed up to 12° incidence angle. For incidence angles larger than 12°, Bragg scattering produced by 
small-scale waves becomes more significant, and thus, two-scale models have to be used (Barrick 
1974; Brown 1978). Similarly, Fig. 8c shows radar measurements and model results using the 
surface slope relationship from TRMM data. For the same surface wind speed, TRMM’s surface 
slope relationship predicts higher a 0 for low incidence angle (<10°), and more rapid decrease of a 0 
along with an increase of the incidence angle. For this case, radar measurements are closest to the 
6.5 ms" 1 curve. 

On 26 July 2002, the ER-2 flew a similar flight pattern to the Caribbean sea as on 9 July. 
Figure 9 shows a portion of the flight track near two directional turns, with an insert of the full 
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flight pattern. The aircraft flew an almost identical track on the outbound and return trips. Similar 
to Fig. 8, Fig. 10 shows <7 0 versus incidence angle during the turns close to way point F. For the 
outbound trip, the ER-2 flew through E, F, and G at 1721 UTC, and er 0 is shown in Fig. 10 as “+” . 
During the return trip, the ER-2 flew through G, F, and E at 1956 UTC, and the corresponding 
curve is presented by “0” in Fig. 10, Surface wind speed measured by dropsonde was 2.5 ms' 1 at 
1753 UTC (N18.51 0 , W84.48 0 ), and 2.3 ms' 1 at 2017 UTC (N24.0°, W86.15 0 ). The quasi-specular 
models using different surface slope relationships with 2.5 ms' 1 surface wind speed are shown by 
the curves in Fig. 10. It is evident that radar measurements from the outbound turn and return trip 
turn agree very well even though they were two-and-half hours apart in time. At a low incidence 
angle, radar measurements match well with model results and are in closer agreement with Wu’s 
(1972, 1990) relationship than the other two. It is worth noting that although the models are “tuned” 
using C e , which is estimated from CRS measurements, the shape of a 0 versus 6 curve will not be 
affected by possible errors in C e since these errors only result in offsets in the a 0 . 

Other turn of events from different days of the experiment were also used to calculate a 0 as 
a function of incidence angle. Figure 1 1 shows cr 0 versus incidence angle from 12 clear weather 
cases. Each case is indicated by a different letter. The different dependencies of a 0 with incidence 
angle are due to the difference in ocean surface structure attributed to different surface winds. 
Results from the quasispecular model with three different surface slope relationships and with 1.25 
ms' 1 and 8.0 ms' 1 wind speeds, are plotted for comparison. It is noteworthy that near the 10° 
incidence angle, a 0 measurements are less sensitive to wind speed, which’ also agrees with the 
quasi-specular model and measurements obtained at lower microwave frequencies (Jones et al. 
1976; Plant 1977; Masuko et al. 1986). As discussed in Section 4c, a 0 at 10° incidence angle were 
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obtained from turn events, during which surface winds were within 3.0-10.0 ms' 1 . The mean value 
of the measured cr 0 is 5.85 dB with a standard deviation of 0.6 dB. 

b. Wind speed and effect of direction on a 0 

Figure 12 presents measured o 0 at 3°, 10°, and 15° incidence angles versus ocean surface 
wind speed. Results using the quasispecular model and three different surface slope relationships 
are also shown for comparison. The discrepancy between the radar measurements and model results 
is likely due to: a) the uncertainty of the wind speed estimate caused by a collocation error between 
the radar beam footprint on the ocean surface and the dropsonde fall position at the surface; b) an 
error in the estimates of the water vapor and oxygen attenuation; c) possible errors in the model. 

Measurements at microwave frequencies have shown that not only the wind speed, but also 
the wind direction and polarization of the radar beam affect scattering from the ocean surface. 
Jones et al. (1976) showed that o 0 were slightly different for upwind, downwind and crosswind 
conditions at a low incidence angle. However, this difference becomes more significant when 
incidence angles are larger than 15°. During CRYSTAL-FACE, the CRS transmit polarization was 
such that the E-field was perpendicular to the direction of flight, thereby, the radar signal was 
vertically polarized (V-plane) relative to the ocean surface during the turns. 

For the 26 July case shown above, the ER-2 made nearly a right angle turn at 1700 UTC 
during the outbound flight (through A, B, C and D in Fig. 9). A dropsonde, launched at 2117 UTC 
on the return trip, measured a surface wind speed of 2.3 ms' 1 from 211°. Figure 13 presents o Q 
versus incidence angle during the turn. The radar beam was oriented upward and slightly to the 
right of the wind in the first half turn (shown by “+” in Fig. 13) and then changed over to the left of 
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the wind in the second half of the turn (shown by “0” in Fig. 13). For an incidence angle larger 
than 12 °, the discrepancy of a 0 between the first half and second half of the turn is evident. 


From 1731:38 to 1734:29 UTC on 29 July 2002, the ER-2 performed a 250° clockwise turn 
in azimuth while maintaining a constant roll (-25°) and incidence angle (-29°). Figure 15 shows o 0 
versus the ER-2 heading angle from this turn. The nearest available buoy (FWYF1 at N25.59 0 
W80.10 0 ) measured a 3.3 ms' 1 wind 57° from the north. The radar beam scanned in the cross-track 
plane which was 90° from the ER-2 heading. Therefore, the radar beam was pointing approximately 
upwind when the ER-2 heading angle was 227°. As evident in Fig. 14, a 0 reached its maximum 
when the radar beam was close to upwind and a minimum when the radar beam was pointing in the 
crosswind direction. These observations are in general agreement with the measurements made by 
Jones et al. (1976) at 13.9 GHz. Note that the offset between the cr 0 maximum and the upwind 
direction is likely due to a small wind-direction error since the buoy was about 22 km away from 
the center of the turn. 

c. Comparison of ocean and land surface return 

The CRS-measured surface returns were highly variable over different surface conditions 
(topography, vegetation, etc.). Compared to terrestrial topography, the ocean surface is more 
homogeneous, and thus the ocean surface return is less variable. Figure 15 shows an ER-2 flight leg 
(2022:49 to 2122:25 UTC) on 13 July 2002 covering ocean and land backgrounds, while Fig. 16 
shows the measured sea surface er 0 between points A and C of this line (2056 to 2112 UTC). 
Between points A (2056:46 UTC) and point B (2104:53 UTC ), where the ER-2 was flying over the 
ocean, the standard deviation of cr 0 is about 0.307 dB; The standard deviation for the inland flight 
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portion (point B to point C) is 3.12 dB. This significant increase in the standard deviation of cr 0 is 
due to the irregularity of the topography and possibly speckle reflection from natural or man- made 
objects. Since we know that the radar beam incidence angle did not vary significantly during the 
flight path, the transition point between ocean and land can be found by examining the standard 
deviation of a 0 measurements. By combining this information with navigational data, the radar 
beam pointing angle in the along-track plane at the ocean/land crossing point then can be estimated. 
One of the practical applications of this method is to determine the approximate pointing direction 
of the beam for a spacebome cloud radar. 

5. Conclusion 

During the CRYSTAL-FACE experiment, ocean surface backscattering measurements were 
obtained using the 94-GHz CRS airborne cloud radar. The calibration uncertainty of the CRS is 
estimated to be 1 dB and was derived from measurements using a trihedral comer reflector and by 
intercomparing measurements with another well-calibrated radar (Li et al. 2003). The calibration 
was applied to deriving estimates of the normalized ocean surface backscatter cross section a 0 from 
the surface return. Attenuation due to water-vapor and oxygen absorption is estimated and 
corrected using a Liebe millimeter-wave propagation model and the meteorological data collected 
by the ER-2 dropsondes. The cr 0 measured by the CRS is then compared with a quasi-specular 
ocean surface scattering model using different surface slope relationships. The analysis results 
show good agreement between the theoretical models and the measurements. It was confirmed that 
even at 94 GHz, the quasi-specular ocean surface scattering model works well for a near-nadir 
operational mode. The CRS measurements also confirm that the o a is insensitive to surface wind 
conditions near a 10° incidence angle, a finding similar to what has been found in the literature for 
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lower frequencies. In addition, the CRS measurements show that the dependence of cr 0 on surface 
wind direction is in general agreement with the measurements made at lower microwave 
frequencies. However, more measurements at 94 GHz are necessary to study the details of this 
dependence. Moreover, the comparison between the ocean and inland surface returns suggests that 
the position of the coastline can be accurately determined, which, in turn, can be used for resolving 
the pointing direction of a spacebome cloud radar. The analyses in this paper support the 
proposition of using the ocean surface as a calibration reference for airborne and spacebome 
millimeter-wave cloud radars. 

When applying a quasi-specular model to the calibration of airborne or spacebome cloud 
radars, errors may result from the uncertainties in the water-vapor and oxygen attenuation 
correction, surface wind measurements, and sampling of the surface return. In addition, the 
effective Fresnel reflection coefficient T e (0,X), instead of F( 0,X), has to be used to account for the 
effects of surface roughness. Since the CRS was calibrated and the calibration results were verified 
using different approaches (Li et al. 2003), this enable us to estimate the correction factor C e using 
a 0 measured by the CRS during the CRYSTAL-FACE. In our study, C e is estimated as 0.88 with an 
uncertainty of 0.16, which corresponds to 1.6 dB in the a 0 calculation (1.0 dB due to the radar 
calibration, and 0.6 dB in a 0 measurements). Using this C e and the quasi-specular model, a 2.2 dB 
calibration uncertainty is achievable for a spacebome radar using its a 0 measurements at a 10° 
incidence angle (a 1.6 dB uncertainty in the C e estimate and also assuming 0.6 dB uncertainty in a 0 
measurements from the spacebome radar). More experimental data, such as measurements of 
backscatter from the ocean surface or wave tank using well-calibrated 94 GHz radars, is necessary 
to verify the value of C e for different geographic locations and surface wind conditions, and to 
reduce the uncertainty in the C e estimate. It is worth noting that even though C e is used by the 
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quasi-specular model to estimate er 0 , the shape of a 0 versus 0 curve will not be affected by possible 
errors in C e since these errors result only in offsets in cr 0 . 
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List of Figures 


Figure 1: CRS inter-comparison with CPRS W-band at the University of Massachusetts Amherst, 
December 13, 2004. The mean difference between two profiles is within 1 dB. 

Figure 2. Ocean surface a 0 versus incidence angle and surface wind speed predicted by the quasi 
specular scattering model 

Figure 3. Geometry of the CRS measurements from NASA ER-2 during the CRYSTAL-FACE 
experiment. The CRS was installed in a nadir-looking mode. The radar beam was 
scanned in a cross-track direction when ER-2 rolled its body to one side during turns. 

Figure 4. Dropsonde measurements at 19:54 UTC, 9 July, 2002. (a) relative humidity, (b) 
temperature, (c) pressure. The two-way cumulative path attenuation due to water-vapor 
and oxygen absorption at 94 GHz is shown in (d). 

Figure 5. (a) Two-way path-integrated attenuation due to water-vapor and oxygen absorption at 94 
GHz derived from dropsonde data collected from CRYSTAL-FACE. It shows the PIAs 
with a mean of 5.8 dB and a standard deviation of 0.65 dB. (b) Ocean surface wind 
speed measured by dropsondes during CRYSTAL-FACE. 

Figure 6. ER-2 flight track on 9 July 2002. It shows three dropsonde launching positions near the 
ER-2 turning point B. 

Figure 7. ER-2 navigational data, radar beam incidence angle and ocean surface return power 
during the turn through point A, B, and C shown in Fig. 6. 
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Figure 8. ct 0 measured by the CRS versus incidence angle during the turn shown in Fig. 6 and Fig. 

7. and compared to qusi-specular model using (a) Cox and Munk’s (1952), (b) Wu’s, 
(1972, 1990) and (c) TRMM’s surface slope relationship. 

Figure 9. ER-2 flight track on 26 July 2002. A larger scope plot is shown at the upper left-hand 
comer. ER-2 flew the track on its outbound trip and on its return trip. It made turns at 
way points B and F. 

Figure 10. a 0 measured by the CRS versus incidence angle from the turns near way point F shown 
in Fig. 8. Cross represents for measurements from the outbound turn through E, F, and 
G consequentially. Diamond represents measurements from the turn through G, F, and E 
on the return trip consequentially. 

Figure 1 1. a 0 measured by the CRS versus incidence angle from different days. A total of 12 turns 
from clear weather is shown. 

Figure 12. a Q versus ocean surface wind speed at different incidence angle. Results from the quasi- 
specular scattering model are shown as the solid, dotted and dashed lines. 

Figure 13. c Q from a turn with 100° heading change. It indicates that a 0 varies relative to wind 
direction for incidence angle larger than 12°. 

Figure 14. cr 0 versus ER-2 heading angle. 

Figure 15. ER-2 flight track on 13 July 2002. It crossed the coastline at point B (21:04:53 UTC). 

Figure 16. a 0 measured by the CRS from 20:56 to 21:12 UTC. From the ocean to inland, it shows 
significant increase in the standard deviation of a 0 . 
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Table 1: CRS system specifications during CRYSTAL-FACE 


Frequency (GHz) 

Peak power (kW) 

PRF (kHz) 

Pulse width (ps) 
Transmit polarization 
Receive polarization 
Antenna beam width ( 0 ) 
Antenna gain (dB) 
Sensitivity (dBZe) * 


94.155 

1.7 

4/5 (Dual PRF) 

1.0 

H 

VandH 

0.6x0.8 

46.4 

-29 (from flight data). 


* At 10-km range, 150-m range resolution, 1-sec averaging 
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g 1- CRS inter-comparison with CPrq w k 
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Fig. 2: Ocean surface a 0 versus incidence angle and surface wind speed predicted by the 
specular scattering model. 


quasi- 
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track direction when ER-2 rolled its body to one side during turns. 
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Fig. 4: Dropsonde measurements 
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Day in July, 2002 



Day in July, 2002 


Fig. 5: (a) Two-way path-integrated attenuation due to water vapor and oxygen absotption at 94 
GHz derived from dropsonde data collected from CRYSTAL-FACE. It shows the PIAs with a 
mean of 5.8 dB and a standard deviation of 0.55 dOB. (b) Ocean surface wind speed measured by 
dropsondes during CRYSTAL-FACE. 
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Fig. 6: ER-2 flight track on 9 July 2002. Three dropsonde launch positions near the 
ER-2 turning point (B) are indicated. 
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Fig. 7: ER-2 navigational data, radar beam incidence angle and ocean surface return power during 
the turn through point A, B, and C as shown in Fig. 6. 
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ER2 Flight Track 20020726 



Fig. 9: ER-2 flight track on July 26, 2002. A larger scope plot is shown at the upper left-hand 
comer. ER-2 flew the track on its outbound trip and on its return trip. It made turns at 1 
points B and F. 



CRS <j 0 (20020726) 



Fig. 10: (T 0 measured by CRS versus incidence angle from the turns near way point F shown in 
Fig. 9. Cross stands for measurements from the outbound turn through E, F, G consequentially. 
Diamond represents measurements from the turn through G, F, E on the return trip consequentially. 


39 


CRS (j 0 Versus Incidence Angle 

20 


^ 10 

o 

b 

w 

u 

“ 0 
T3 
<D 
N 

D 

E 

z -10 


-20 

0 5 10 15 20 25 

Incidence Angle (Degree) 

Fig. 11: cr 0 measured by CRS versus incidence angle from different days. Twelve turns from 
clear weather are shown. 
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c T 0 Versus Wind Speed at 94 GHz 



Fig. 12: a 0 versus ocean surface wind speed at different incidence angles. Data points are given by 
symbols. Results from the quasi-specular scattering model are shown as the solid, dotted and 
dashed lines. 
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Fig. 13: < 7 0 from a turn with 100° heading change. It indicates that <r 0 varies relative to wind 
direction for incidence angle larger than 12°. 
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Fig. 14: o 0 versus ER-2 heading angle. 
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Fig. 16: <r 0 measured by CRS from 2056 to 2112 UTC. From the ocean to inland, it shows 
significant increase in standard deviation of <x 0 . 
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